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ABSTRACT 


Simultaneous measurements of mean wind speed, horizontal wind 
velocity, vertical wind velocity, and water height were made during the 
summer of 1968 at an exposed field site off the New England Coast. A 
vertical array of five cup anemometers provided wind speed data for 
heights extending up to 8.2 meters. Hot film anemometers operated 
in pairs measured component wind velocities at a height of 1 meter 
above the mean water level. Hot film anemometers mounted on two 
horizontally separated instrument platforms measured a two-point 
spatial correlation field. Measurements were conducted with and with- 
out an artificial sea slick on the water, demonstrating the importance 
of small waves to air-sea interaction. 


Analysis of the wind profile data indicates better than 90 percent 
of the validly measured wind profiles are logarithmic. The friction 
velocity Ux is a linear function of wind speed in the range of wind 
speeds studied, while the roughness length parameter, z,., is highly 
Structured. The structure of the roughness length parameter can be 
explained in part by the Kelvin-Helmholtz instability and the onset 
of wave generation. A discoutinuity at 4 meters/sec~+ cannot be a 
explained. The drag coefficient is constant with a value of 1.6 x 10 
except in the vicinity of roughness length discontinuities. The mean 
wind speed at 10 meters is represented as the sum of a mean profile 
term and a velocity defect term constructed from the correlation 
between Ux and z,.. The velocity defect term has a systematic behavior 
as a function of wind speed. 
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Taylor's Hypothesis is tested based on cross-correlations and 
cross spectral denisty computations of the data from spatially separated 
hot film anemometers. Taylor's Hypothesis does not appear to be consis- 
tently valid for atmospheric turbulence in the general spectral domain 


of ocean surface waves. Reynolds stress measurements are consistent 
with the values cbserved by other investigators and with the values 
inferred from the wind profile data. Cross spectral density computa- 
tions indicate that the near-surface horizontal velocity is a maximum 
over the wave trough for both the well and the wind waves, except 
when the water is covered by a slick. In the case of a slick, the 
horizontal velocity has its maximum value over the wave crest of the 
swell and looses coherence with the wind sea. 


Thesis Supervisor: Erik L. Mollo-Christensen 


Title Professor of Meteorology 
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Introduction 


This paper is concerned with the momentum transfer processes at 
the air-sea interface; specifically those processes which contribute 
to the growth and maintenance of wind waves. Attention is focused on 
the atmosphere over the ocean, and the isolation and identification of 
atmospheric phenomena unique to the ocean atmosphere. As a matter of 
viewpoint, the research was approached with a minimum concern for the 
air's action on the sea, as opposed to the sea's action upon the atmos- 
phere, in the belief that the atmosphere's ability to transfer energy 
into the ocean is dependent primarily on the nature of the effects of 
the ocean on the atmosphere. 

As a point of departure, we look at the atmosphere directly 
over the ocean from the point of view of the atmospheric investigator 
as contrasted with that view commonly held by the wind wave investigator. 
The former tends to consider the ocean surface as just another rigid 
boundary, where the magnitude and direction of the energy fluxes across 
the boundary are not considered. [Paulson (1967); Pond, Smith, Hamblin, 
and Burling (1967); Pond, Stewart, and Burling (1963).] While such an 
approach appears to result in the proper wind profiles over the sea, 
it is obviously the incorrect viewpoint when discussing the energetics 
of the air-sea interface. Stewart (1967) has shown that the saturated 
wind wave field has gained momentum from the atmosphere equivalent to 
the entire momentum contained in the atmospheric mean wind extending 


to a height of one ocean wave wavelength (roughly 40-60 meters as a typical 








figure). Intuitively, during a falling sea, one should expect the 
momentum flux across the air-sea boundary to reverse sign. 

The wind wave investigator, looking up from below, tends to treat 
the atmosphere as an infinite energy source over the ocean, being milked 
of momentum in response to proper excitation from the wave below. Again 
this appears to be a correct viewpoint in a most general sense, however 
one is left asking just how much momentum can the atmosphere provide. 

The focal point for the different viewpoints in the treatment of 
the atmosphere is the vertical mean wind profile. Stewart (1961) and 
Kinsman (1965) argue that the magnitude of the momentum flux from the 
atmosphere to the ocean is so large that the atmosphere should not be 
able to support a logarithmic profile. Atmospheric investigations, 
such as those summarized by Roll (1964), show that the oceanic profile 
is in general a logarithmic profile, but beg the question of explaining 
the maintenance of the profile in the face of wave generation. All of 
this suggests to this investigator that the oceanic vertical mean wind 
profile must be considered as the result of a quite different process 
from that effecting its rigid boundary counterpart. 

The research reported in this paper has tried to identify unique 
aspects of the oceanic wind profile, and there find some explanation 
for the existence of these unique aspects. The research, and consequently 
this report, can be conveniently divided into three general areas. The 
first area is the development and use of suitable measurement techniques 
needed to observe the atmosphereover the ocean. Next, the vertical 


mean wind profile is examined. Finally, some information is presented 
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concerning the two-dimensional motions directly above the ocean. As an 
illuminating side bearing on the results, some experimental results 


with artificial sea-slicks are reported. 
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Site and Instrumentation 


Field measurements were conducted from a rigid spar designed and 
erected under the supervision of Professor Erik Mollo-Christensen in 
Buzzard's Bay, Massachusetts. The spar was placed in 60 feet of water 
as shown in figure 1.1, and provided an ideal platform for mounting 
sensitive instruments for air-sea interaction studies. 

Figure 1.2 is a rough sketch of the geography of the area around 
the spar. Prevailing winds for the area are from the Southwest, and 
almost all of the observations recorded were of winds from the Southern 
sector. The nearest land to the Southwest is Block Island, at a distance 
of 26 miles. The bottom topography to the South and Southwest gradually 
increases in depth, so one need not be concerned about wave field distor- 
tion due to upwind and up-sea shoaling of the wind-waves studied in this 
report. 

The deep water, the high wind speed variability, the virtually un- 
limited fetches to the South, and the rigid instrument platform provide 
an excellent environment for the study of air-sea interaction processes. 

Figure 1.3 is a photograph of the spar with instruments rigged. 
Note that the instruments are mounted on the end of a trainable instru- 
ment boom extending 12 feet from the spar. The boom was oriented into 
the wind at the beginning of each run, and was elevated or depressed 
so that an inscribed boom "waterline" was positioned at the mean water 
level of the sea expected for the duration of the data run. During any 
given data run, the mean water height would vary less than 1] foot due to 
tidal effects, and was treated as being essentially constant during the 


run. 
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Figure 1.1 


Bottom mounted spar design for use in 
Buzzards Bay during the summer of 1968. 
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For spatial wind correlation studies, a taut moored, moveable, wind 
stabilized instrument buoy as placed at various locations near the main 
spar. Initial plans called for two such buoys, as shown in figure 1.4, 
however, one of the buoys was lost during heavy weather before experiments 
could be conducted using two buoys simultaneously. The portable buoy 
consisted of an instrument cage mounted on the end of a 25 foot section 
of 2 inch structual aluminum pipe, of which 6 to 10 feet protruded 
above the water. The floatation for the buoy was roughly 10 feet below 
the water line. The instrument electronics were contained in a Sears- 
Roebuck pressure cooker, as shown in figure 1.5. A hot film probe was 
mounted on the end of a two meter long combination instrument support 
and radio antenna. The entire instrument package was dropped into bearing 
sleeves mounted in the instrument cage when the buoy was to be used, 
otherwise the instruments were maintained ashore. 

The instrumentation on the spar and the portable buoy remained 
unchanged during the entire summer of 1968, except for minor variations 
to accommodate special experiments. The instrumentation consisted of: 

One Beckman and Whitley 6-cup anemometer and wind vane 

System mounted at the top of the spar at a height of 8.2 meters. 

Four C. W. Thorthwaite Associates cup anemometers located 

at 1, 2, 3, and 5 meters above the mean water level. 

Two sets of DISA S & B battery operated hot film anemometers 
mounted at 1 and 5 meters above the mean water level. The hot 

film probes were used in pairs as described in chapter II, to 


provide direct measurement of Reynold's stress. 








Figure 1.3 Instrument spar with boom ex- 
tended and instruments rigged. 








Figure 1.4 Satellite Buoy. 
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One capacitance wave gauge designed by the Department of Meteor- 
ology Fluid Dynamics Laboratory, mounted directly below the hot film 
anemometers. 

One set of slope-matched thermistors mounted at 1 meter and at 
5 meters to give temperature difference readings for the computation 
of Richardson's number. 

One hot film anemometer and associated circuitry located on 
the portable buoy. 

The signals from the various sensors on the spar were preprocessed 
as necessary, and then transmitted through a floating umbilical cord 
to the R. R. Shrock, moored by a two point moor 150 feet downwind 
from the spar. The signals from the portable buoy were transmitted 
by radio link. The data logging system on board the Shrock is shown 
in figure 1.6. The Thornthwaite Cup anemometer data was recorded 
and printed out by a Thornthwaite model 106 wind profile register 
system. The temperature data and the wind direction was recorded 
on strip-chart recorders. The remaining data was recorded on a 
Precision Instrument 6208 8-channel tape recorder. 

Considerable attention was given to the transmission of the 
hotwire signals. The outputs of these anemometers were linearized 
using an in-house designed linearizer (see appendix 3) and then 
converted into a frequency modulated signal prior to transmission 
and recording. Such preprocessing is considered essential to 
maintaning signal fidelity and low signal-to-noise ratios for the 
hot film signals. In the case of hot films operating together as 


pairs for the direct measurement of Reynolds stress, an additional 








Figure 1.5 





Satellite buoy instrument package 
constructed from pressure cooker 
with instruments in place. 
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Figure 1.6 Data logging system on board the R/V R.R. Shrock. 
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Figure 1.7 Signal conditioning amplifier (left) 
and voltage-controlled oscillator (right) 
constructed for use with hotwire systems. 
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precaution was taken by multiplexing the two signals prior to trans- 
mission and recording. This avoided any question of phase match 
errors due to multipath transmission effects, tape recorder head 
alignment, etc. 

With the exception of the temperature difference measurements 
and the wave gauges, all of the equipment consisted of standard 
off-the-shelf hardware. A full description of the cup anemometers, 
the wave gauges, and their characteristics is covered by Seesholtz 
(1968). A discussion of the Reynolds stress measurement is contained 
in Chapter II. 

Turbulent fluctuations in the atmosphere were measured using 
DISA type 55D05 hot wire anemometers in conjunction with DISA 55A85 
hot film probes. The anemometers are battery operated, compact, and 
ideally suited to the rugged field environment. The probes consist 
of a platinum film sensing element fused to a wedge-shaped tip of a 
quartz support. The film is electrically insulated from the salt air 
by a thin quartz coating over the film. The probes were operated at 
an overheat ratio of 1.7 to 1, recommended by the manufacturer as 
being the optimum for a probe operated in air and subjected to dousing 
from occasional rain drops or intermittent water immersion. While the 
probes performed well in the field and survived environmental extremes, 
they were subject to sudden and unanticipated step changes in their 
characteristic operating resistance, requiring close monitoring. 
several days of promising data had to be discarded because of this 
probe behavior. Additional probe characteristics are discussed in 


Chapter II. 
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Figure 1.7 shows the conditioning amplifiers and voltage controlled 
oscillators employed with the hotwire anemometers. The amplifiers were 
built around Texas Instrument integrated operationalamplifier chips. 

The voltage controlled oscillator is a modification to the standard 
design presented in Navships 93484 (selected semiconductor circuits). 

The temperature difference measuring device consisted of two rug- 
gedized bead thermistors encased in glass, each mounted in identical, 
naturally ventilated radiation shields. The thermistor beads were slope- 
matched to provide an identical resistance variation for a given temper- 
ature variation to within one part in 10>. Each thermistor was connected 
as opposing elements of a balanced bridge, the output of which was 
amplified by standard operational amplifiers and calibrated to provide 
a full scale meter deflection on board the R. R. Shrock chart recorder 
in response to a 1 degree centigrade temperature difference. The high 
mass of the glass around each bead resulted in a system time constant 
of approximately 20 seconds, providing high damping and stability. 

After the system had settled down, it was capable of measuring temper- 
ature differences to with .01 degrees centigrade in an ideal environment. 

While the instrument worked quite satisfactorily in the laboratory 
and during pierside tests, it was somewhat of a disappointment at sea. 

As shown by the tabulated Richardson number calculations in Appendix I, 

the instrument would work well when the sea conditions closely approximated 
the ideal conditions available at pierside. Otherwise it would not work 

at all. The primary problem was the wet bulb effect caused by the depo- 
sition of salt spray on the bead. A suitable solution was not found for 


this problem. 
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Chapter II 


The Direct Measurement of Reynolds Stress 


General. 

The direct measurement of Reynolds stress at some point near the 
water surface has been one of the major problems in the observation of 
air-sea interaction processes. Instruments capable of providing such 
measurements are usually quite delicate, and require special handling 
and care to survive the hostile marine environment. The nature of 
the physical processes under study further impose demanding require- 
ments of dynamic response and sensitivity far beyond that required of 
the laboratory test environment. The importance of Reynolds stress 
measurements to the understanding of the near ocean atmospheric tur- 
bulence, coupled with the lack of literature discussing the instrumental 
problem, requires that this problem be separately considered in detail. 

The over-riding instrumental requirement in air-sea interaction 
studies is resolution, with dynamic range running a close second. Wind 
perturbations induced by the underlying wave action represents a small 
percentage of the overall turbulent energy. Consequently the instruments 
must not only have the dynamic range to handle the large gusts charac- 
teristic of the atmosphere in general, but must also have the sensitivity 
to detect the small wave-induced fluctuations. Hot film probes, operated 
in orthogonally mounted pairs, have been chosen as being best suited for 
this study. They fulfill the requirements of sensitivity and dynamic 
range, and have the rugged construction essential to instruments oper- 


ated in the maritime environment. 
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Instrument Calibration and Response 


From the outset it is assumed that the Reynolds stress sensor is 
mounted on a rigid mount near the ocean surface. Pond (1968) has dis- 
cussed the effects of buoy motion, which must be considered if a rigid 
mount is not available. The problem of Reynolds stress measurement in 
general will be treated first, followed by the specific case of a hot 
film probe pair. 

Kraus (1968) has shown for an instrument whose alignment deviates 
from a true vertical by an angle of one degree that one can expect 
errors in the measurement of Reynolds stress in excess of 100 percent. 
Kraus suggested several indirect approaches which might be suitable for 
correcting this error, both of which are considered unsuited for near- 
ocean measurements. An alternate approach to those of Kraus is to 
invoke the continuity condition that the mean vertical velocity near 
the boundary goes to zero, and use this to correct for instrument align- 
ment error. 

Let U and W represent horizontal and vertical components of the 
wind. The apparent velocities as sensed by the instrument will be 
denoted by a subscript (a). Following the usual Reynolds convention, 
an overbar denotes a mean quantity, while a prime represents a deviation 
from the mean. From the geometry of figure 2.1, one can write: 

Wa + U sina + W cosa 
Ua = U cosa —- W sina (2.1) 
Applying the Reynolds convention and invoking the requirement that 


W = QO, equation (2.1) can be written: 








Zsa 





Z 
fa 
\ 
\ 
\ 
\ 
\ 
\ 
\ 
\ 
\ 
\ 
\ 
\ 
\ 
\ 
\ 
\ 
\ 
\ 
‘ 
\ ae Xa 
\ ee 
\ _-7 INSTRUMENT 
\ Be AXIS 


x 


TRUE AXIS 


FIGURE 2.1 


Geometery of Instrument alignment error. 
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Wa = U sin QO 
Ua = U cos @ 
t t 


W a= WU %sin® + w ocasa 


u cosa - w' sina (2.33 


& 
a) 
II 


Forming the Reynolds stress term from equation (2.2): 








u'w' = U'aW'a cos 2a - Ssin 20400h 4 = Wet ) 
(233) 
Equation (2.3) can be completely evaluated provided the angle a can be 


determined. From equation (2.2), however, we can write: 


(25G)) 


Sn Ge mam 
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One can use equation (2.3) and equation (2.4) together to compute the 
Reynold stress, since one can directly measure the apparent vertical 
velocity and the mean horizontal velocity. If the two equations are 
combined into a single expression, and the small angle approximation 


, 2 
is made to order Q,;: 


2 = 2 2 
Wa(U'. | we ) 


ewe = UAW [1-2 ] 
U U 5) 





For processing of hotwire signals, however, it is more convenient 
to express the equations directly in terms of the hotwire output signals, 
rather than the derived expressions. An ideal hotwire would have a 
cosine law response, where the output signal from the hotwire varies 
directly as the cosine of the angle between the wind direction and a 
normal to the hotwire. In reality, however, the hotwire response 
deviates from the ideal cosinse law response due to the edge cooling 


of the wire. A better approximation for the directional response of 
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a hotwire is modified sine-cosine law response of the form: 


E, = (Vcos9+KV sin 9) (256) 

E, is the voltage output of the hotwire aenometer, and at this point 
will be assumed to be a linear function of flow velocity. The non-lin- 
earity of the hotwire set putput will be treated separately. V is the 
free stream velocity, and angle 9 is the angle between a normal to the 
hotwire and the free stream velocity vector. K is an experimentally 
determined constant characterizing the hotwire response. For the hot 
film probes used to measure Reynolds stress, K = Q.55. 

Consider two hot film probes mounted as shown in figure 2.2. The 
geometry of the probes is shown in 2.3. From this geometry and from 


equation (2.6) one can write: 


Vp = V (cos (45 + 0) + K sin (45 + 0)) 


Vo> V (cos (45 — ©) + Kk sino. ae) (2.7) 

Vp and V, represents the components of the free stream velocity sensed 
by hotwires oriented with the wires each normal to the p and s axis, 
respectively. Also, from figure 2.3: 


U Vees 9 


W 


Equation (2.7) can be manipulated by expanding the trigonometric functions, 
and then forming new expressions from the sum and difference of equation 


(2.7). After evaluating known angles, the result becomes: 
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Sketch of hot film probe arrangement. 





FIGURE 2.35 


Geometery of hot film probe arrangement. 
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FIGURE 2.4 


Orthogonally mounted hot-film probe directional 
response based on wind tunnel tests. 
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The response characteristics of a hotwire pair mounted in an 
orthogonal configuration is shown in figure 2.4. Note that U and W 
components are accurately measured within plus or minus 30 degrees of 
tue horizontal. This requires that lwl<o.6|ul, a condition which is 
generally true in the atmosphere since W = 0 and w* is at most an 
order of magnitude smaller than U. 

The probes encompass a field of measurement of roughly one-half 
centimeter, while the scale of the fluctuations being studied are of 
order 10 centimeters or greater. Consequently the measurement field 
of the probes can be considered to be a point. Further, for measurement 
of mean wind speeds over time intervals of several minutes, one can 
expect two sensors separated horizontally by a distance of roughly l 


meter to sense the same mean wind speed. Therefore: 


oe Ve 
yz «yz 
Pp S 
i i ve 7 oie (2.10) 


Where Vs is the horizontal mean wind cwveed as sensed by a cup anemometer 


located at the same height and in close proximity to the hotvvire pair. 


B is an arbitrary constant of proportionality to be determined. Equation 


(2.10) says that both probes and the 


+ 
t 
ft 


up anemometer sense the same mean 
wind, while each of the two probes sense the same mean scuare turbulent 
fiuctuations. Reconsidering the alignment problem, if the output of 
the two probes is adjusted, to satisfy equation (2.10), then W, = 0 and 
the error terms in equaticn (2.3) are eliminated. 

Given two records of hotwire data from a pair of hotwire sensors 
mounted for Reynolds stress measurement, the records can be statistically 


normalized with a mean equal to zero and with a unit standard deviation. 
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This satisfies equations (2.10) and (2.3). The signals would have an 
mtoaerary calibration, however, and the velocities could not be deter- 
mined arilees | is determined. To determine B we will depend on anneen 
wind trend within the wind record. Let 

MV. = e 

MV, ~ Ue (2) 
where M is a velocity calibration constant; i and Vi. are the values 
of velocity at any given instant in the two sets of normalized velocity 
records. From equation (2.10) and equation (2.11) one can write an 
expression for a record possessing a mean velocity trend. The mean 


velocity over an arbitrary sub-interval within the record about times 


ty and t, can be expressed as: 


M(V_(t,) te V(t) 


BV (ty) 

M(V_(t,) af V(t 2)) = BV (t.) CZ alee) 
Where the two time intervals are non-overlapping, and by virtue of the 
trend, Vi (ty) is not equal to. V (to). Adding the two expressions in 


equation (2.12) and solving for B/M: 


B_1 _ (Va(t) - VaCty)] + [Vpep) - Weep) 
‘i [V.(tp) - Ve(ty)] 


co 


This expression can be used to scale U by substituting into equation 


(2.9) with the result: 
= Vv. 14 
ah Ve + V,) Es . (2.14) 
The scaling of the vertical wind component is based on an examinacion 


of the followirg equations, which follow directly from equations (2.9) 


W7eloy and (2.14): 
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If we multiply out both terms, and substitute an expression for U'” 


where appearing in W'* , we can write: 


2 2 
wees aye Ww - Vv)? (2.15) 
1 


Finally, expressing Ny and No in terms of K, the result becomes: 


er CK 


W=AG—) W, - ¥) (2.16) 


This equation, along with (2.14), has been used to solve for U and W 
using hotwire signals. These equations incorporate corrections for instru- 


ment alignment, in-site caiibration, and probe response characteristics. 


Other instrument errors. 

There are three classes of errors which will be discussed next. 
The first concern, of course, is the non-linear response characteristics 
of the notwire anemometer. Tne usSuai approach, which was followed during 
the summer of 1968, is to employ an analog device such as the one described 
in appendix III to linearize the anemometer output. In retrospect, how- 


ever, this appraoch was a poor one since it tended to maximize another 


qj - is . eh a, se = : : aon ~~ 7? er ae 
CuiaSS Cr error sz COnRCSETN, i.2., Cne Gala CLaGRSMLSSiOn aia Conversion 


“ 


Cieor. 
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The initial philosophy, based in part on the then existing data 
analysis techniques, was to linearize the signals in the field before 
transmission and recording. The justification for doing this was to 
record and process signals which could be subsequently treated as linear 
signals by analog data processing systems. As the techniques for hot- 
wire data processing were refined, coupled with a shift from analog to 
digital processing, it became clear that the problem of linearization 
needed reexamination. 

The basic problem is that while most hotwire anemometers have 
essentially unlimited dynamic range when used for atmospheric measure- 
ments, the data recording and conversion devices have limited dynamic 
range and limited sensitivity. For example, a 10-bit analog-to-digital 
converter was used to digitize the signals. The maximum range of this 
device is plus or minus 512, for a maximum resolution of roughly one 
part in 1000. If one starts allocating the available range of the digitizer 
to accommodate the maximum wind gust expected of say 20 meters-sec , 
then for a linear signal the maximum resolution is of the order of .Q2 
eeeccec This is the approximate magnitude of the velocities we 
wish to resolve. Consequently significant data is degraded or completely 
lost due to transmission range resolution. 

The response of a hotwire is usually expressed: 


: 2 = gyl/2 (2.17) 


Where E is the hotwire voltage output, E, is the no-flow voltage output, 
8 is an instrument characteristic parameter, and U is wind velocity. 


Note that the hotwire voltage output varies as the fourth root of the 
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wind velocity. Consequently for low wind velocities the set sensitivity 
is large and diminishes for higher wind velocities. The characteristic 
response of the hotwire anemometer therefore provides a signal compression. 
If the E, voltage is subtracted from the signal before transmission and 
recording, one can then utilize the entire channel range for the non- 
linear signal. After transmission and conversion, the signal can then 

be linearized in an analog or digital computer, where there is again 
essentially unlimited range and resolution. The result is a signal with 


wide dynamic range and with a non-linear resolution. The signal will 


have high resolution at low wind velocities --- where it is needed --- 
and poorer resolution at high wind velocities -- where it is of less 
importance. 


Finally, a word about the dynamic response of hot film probes. 
Rasmussen (1968) has shown that the DISA 55A81 pyrex-backed hot film 
probe had dynamic response characteristics markedly different from that 
inferred from static calibration. In view of the unusual geometry of 
the 55A83 and 55A85 probes, dynamic calibration runs were conducted 
using the low-turbulence wind tunnel at the Massachusetts Institute of 
Technology. The results are shown in figure 2.5. The probes have a 
resonant point at about 1.5 Hz. The response falls off to about 50 
percent of the static response at frequencies above 50-100 Hz. The 
region of interest for atmospheric measurements is the region below 0.5 Hz. 
In this region, the dynamic characteristics are well approximated by 


the static calibration of the hot film probes. 
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Chapter III 


The Vertical Mean Wind Profile Over the Ocean 


For Light to Moderate Winds 


The nature of the vertical mean wind profile over the ocean remains 
one of the central issues in the study of the atmospheric boundary layer 
over the ocean. Consequently an experimental study has been undertaken 
in an effort to increase our knowledge of processes in this regime through 
observation. Over a period of two years, measurements of mean wind speeds 
at vertically separated stations have been collected by Seesholtz (1968) 
and this author. This effort has resulted in a data set of 299 mean 
wind profiles measured in the first ten meters of the atmosphere. The 
profiles have a distribution as a function of the wind speed at ten 
meters, taken as the reference height wind speed throughout this report, 
as shown in figure 3.1. The results of the study are valid in the light 
to moderate wind speed regime. 

Each profile is based on wind speeds measured simultaneously from 
four to eight anemometer stations. The lowest station was within one 
meter of the mean sea surface unless sea conditions required otherwise, 
while the uppermost station varied in position between 8 and 12 meters 
above the mean sea surface. While instrument configurations varied, in 
every case the anemometer stations were positioned more or less logarith- 
mically with instruments concentrated within the first few meters above 
the water's surface. The instruments were attached to a rigid spar stand- 
ing in 60 feet of water, therefore, the measurements are free of effects 


such as buoy or ship motion while the behavior of the wind-waves can be 
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expected to be that of deep water wind waves. All wind speeds have 
been averaged over a ten minute interval. 
The basis for the analysis and discussion is the logarithmic wind 


profile, defined as: 


Ux 


UG) = &n(—) aa 


Zo is the roughness length, defined as that height at which U(z), the 
mean wind at height z is equal to zero. This parameter is normally 
thought of as characterizing the effect of boundary on the mean air 


flow. U, is the friction velocity, where U, = Catwryl/2, 


K is von 
Karman's constant, taken to be equal to 0.42. From an experimental view- 
point, Ux and z, can be uniquely determined from equation (3.1) provided 
the profile is in fact logarithmic. These parameters become indeterminate, 
however, for experimental data representing non-logarithmic profiles un- 
less one has a prior knowledge of the functional form of the profile. 
Therefore, experimental profiles which are not logarithmic will be dis- 
carded. Non-logarithmic profiles can result from a myriad of causes, 
Such as stability changes, a passing ship upwind of the observation . 
Station, etc. Aside from the indeterminancy problem, discarding anomalous 
profiles excludes extraneous information which might otherwise bias the 
results and their subsequent interpretation. Fortunately, while valid 
cases of non-logarithmic profiles exist over the ocean, they are the 
anomaly. In this study, thirteen percent of the data sample has been 


discarded. Of that 13 percent, better than half of the cases can be 


discarded on valid physical grounds. Therefore, better than 90 percent 
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of the valid measured profiles are logarithmic, as shown in the shaded 
areas of figure 3.1. 


Equation (3.1) is linear in &n(z), and can be written: 


Qnz = 


(2 eee (3.2) 
U, | 


One can use linear regression techniques to find the best least-squares 
fit of the assumed profile form to the observed data. Following Hoel 
(1954, p 127), equation (3.2) can be written in terms of the observed 


wind speed V at instrument height Y;: 


_ V(X) anY, - (enY, int W (Y5)) 
wae "lav, ) = ee en D 
V(x4)? - VX;) 
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The overbar represents an ensemble mean, while N is the number of pro- 
file points measured for each profile. The ensemble mean is defined in 


the usual sense for discrete data: 
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From equation (3.3) the equation of the line best fitting the profile 
» 
is determinec. From this line, Ux and z, can be directly computed, 


along with a linear correlation coefficient: 
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VY, ) &n¥ = (2nY,, ) (V(X, )) 
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[ (2ny , )2- (&nY,, ) ] Wes) - V(Y,) J 


(3.4) 
These correlation coefficients measure how well the logarithmic profile 
form represents the observed profile. Detailed information for each 
profile is presented in Appendix I. Those examining this appendix will 
note that Ux has been computed for several heights, and does not always 
appear as a constant with height. This variation is because the Ux(z) 
displayed in appendix I is computed from equation (3.2) using the obser- 
ved wind at height z, instead of the computed profile wind. The ten 
meter reference U,, however, in all cases is the profile value. 

Table 3.1 lists the distribution of correlation coefficients com- 
puted from each of the observed mean wind profiles. A correlation coef- 
ficient of 0.94 has been chosen as the division point between those pro- 
files defined as logarithmic and those profiles discarded as deviating 
from the logarithmic form. The shaded region of figure 3.1 represents 
the distribution of the logarithmic wind profiles after the anomalous 
profiles had been discarded. Note that those profiles taken to be log- 
arithmic are distributed as a function of 10 meter wind speed roughly 
the same as the distribution of the total sample. Each of the discarded 
profiles were individually examined, with the following results: 

A. Of the four profiles with correlation coefficients less than 


0.74, three were measured during a time period when the mean wind speed 


suddenly increased at the upper levels during the averaging period. 
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The fourth profile could not be classified due to insufficient supple- 
mentary information. 

B. The profiles with correlation coefficients between 0.94 and 
0.80 consists of 19 profiles measured during experiments with artificial 
sea-slicks. These are not considered representative of a true marine 
environment. The remaining cases, though not attributed to any identi- 
fiable cause, were charatterized by intermittency. That is to say for 
a set of consecutive profiles, there would be one or two anomalous pro- 
files dispersed among an otherwise fully logarithmic set. 

After the profiles were computed and sorted those profiles defined 
as logarithmic were grouped into classes according to ten-meter wind 
speed. While the data was computed for a number of class intervals, a 
class interval of BOC eeu was chosen as optimum for displaying the 
data. A broader class interval retains the gross features, but tends 
to smear the detail. Narrower class intervals result in a discontinous 
histogram. The class boundaries fall on the integral wind speed value 
and the one-half peeer one values of reference wind speed. An obser- 
ved wind speed which fell exactly on the boundary was grouped with 
the higher adjacent class. The data in each class was used to compute 
a mean friction velocity and roughness length appropriate for the 
class. 

The roughness length as a function of reference mean wind speed 
is shown in figure 3.2. The structure of this parameter was quite 


unanticipated. The first appearance of this structure resulted in extensive 
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Figure 3.2 Roughness length of a function of 10 meter wind speed. 
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recomputation, however the features persisted. Wu (1969) recently has 
published a paper with data processed in a somewhat similar fashion. 
While Wu did not take note of it, his roughness length plot shows peaks 
Similar to those in figure 3.2 at 4 and 8 Tee een Wu's graphs 
were drawn without any data discrimination and for wider class inter- 
vals, therefore, the peaks did not stand out as markedly as they do in 
figure 3.2. Additional corroborating evidence of an indirect kind is 
provided by Kinsman's (1968) review of the Beaufort wind scale. The 
scale is shown in figure 3.5. Note that the natural wind speed division 
points selected on the basis of the appearance of the sea surface to 
seamen coincides closely with the peaks in the roughness length graph. 

On the basis that marked changes in roughness length characterize 
a change in the condition of the underlying surface as seen by the at- 
mosphere, one can speculate on the significance of each of the peaks. 
If an air stream was introduced over a perfectly calm sea, and then the 
wind velocity was slowly increased, there should exist some velocity at 
which the first wind-waves will appear. Below this velocity the atmos- 
phere should be expected to see a smooth aerodynamic boundary. Above this 
velocity, the atmosphere should feel a boundary with different aerodynamic 
characteristics. In figure 3.2 such a change occurs at 2 een ee eee 
Suggesting this as the minimum wave generating wind speed. This velocity 
is equivalent to a one meter wind speed of 1.5 Roe eae and is in 
close agreement with Jefferys' minimum wind speed for wind gemration as 
discussed in Lamb (1932, p 625). 

The roughness length peak at 4 meters/ seca cannot be explained. 
To the knowledge of the author, there isn't any theoretical work which 


could be used as a basis for explaining this peak. A tentative physical 
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explanation will be suggested later based on an analysis of the velocity 
defect graph, figure 3.5. The marked peak at 8.5 Woe acre however, 
is reminiscent of the widely discussed Kelvin-Helmholtz shear instability, 
predicted by Munk (1947) to occur at peieeasee = The higher observed 
wind velocity for the onset of the shear instability is attributed to 
the modification of the surface tension of the sea water due to dissolved 
sea salts and contaminants. 

Friction velocity as a function of reference wind speed (figure 3.3) 
yshows a strong linear dependence upon velocity. For the Q-10 meter/sec 
speed range, the friction velocity can be approximated by the formula: 


Uy, = .04U(10) (3.5) 


Note that the friction velocity plot shows only minor deviations from its 
trend in the vicinity of 2 and 4 meter nee 2, At 8.5 neers) saan however, 
the friction velocity appears to be strongly influenced by the same pro- 
cesses which effect the roughness length parameter. This strong effect 

is again what one might expect in the case of the Kelvin-Helmholtz insta- 
bility, where ae boundary instability changes the character of the entire 
turbulent regime. 


lite drag coeificient at height 2 ts defined. 


Ux 2 
C)(z) = Gey (3.6) 


The plot of the drag coefficient is displayed in figure 3.4. Again, 

the singularities are clearly evident, as is the constant drag coeffici- 
ent Obel Oo X ie between singularities. The behavior of the drag 
coefficient in the eae oe the singularities is typical also of the 


behavior of the values of the friction velocity. This wide deviation 
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Figure 3.3 Friction velocity as a function of 10 meter wind speed. 
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Figure 3.4 Drag coefficient as a function of 10 meter wind speed. 
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at the singularities suggests that the variables have a modal structure 
within the class interval. This is what one would expect for sharp dis- 
continuities. Such a structure would result in a mean for the class 
interval lying somewhere between the two modes of the variable, and 
extremely large and somewhat meaningless variance value. The functional 
dependence of the drag coefficient on wind velocity has been a matter 

of extensive debate, based and scattered data. The results presented 
here suggest that the inferred functional dependence of the drag coef- 
ficient would depend more on how fast the wind was blowing during the 
various investigator's measurements or how the data was grouped, than 

on any physical relationship. The value of 1.6 X he for the drag 
coefficient agrees with laboratory values of the drag coefficient deter- 
mined by eddy correlation techniques as summarized by Roll (1965, 

p 158) and the values measured in the field using eddy correlation 
techniques reported by Weiler and Burling (1967). 

The computed mean values of Uy and Z, for each class interval 
were used in equation (3.1) to compute a mean 10 meter wind speed for 
each velocity class. It was expected that the computed speed would 
closely approximate the specified mean wind speed for the class. In 
most cases, however, this computation resulted in a consistent under- 
estimate of the wind speed approomiate for the class. The magnitude 
of this underestimate is shown in figure 3.5. Since we have been pro- 
cessing the data using statistical approaches, this result should be 
interpreted from a statistical viewpoint. The implicit assumption in 
equation (3.1) is that Uy, and Z, are each independent parameters required 


to describe the flow. If, however, they are not statistically independent, 
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the mean profile form can be expected to deviate from that inferred from 
equation (3.1) by the amount of correlation between U, and zo. We will 
represent U, and 2, as the sum of mean quantity and a deviation from the 
mean: 


Ux = Uy + Uy! 


an = = (in = ) (&n >)" (Gan 


° 
Equation (3.7) is substituted into equation (3.1) and an ensemble mean 
is taken. The results are: 
1 


(n>) +7 U,' (in = >! 


= 1 
Uz) = FU K Us = (3.8) 


k ~* 
Note that this is not the same type of averaging process as is done in 
the familiar Reynolds convection. In this case we are averaging over 
an ensemble defined by all those values of U, and z, which fall within 
a specified class interval, rather than over an arbitrary time period. 
Equation (3.8) therefore indicates that for any given class of wind pro- 
file data, the mean wind speed at a height z can be represented as the 
sum of a mean profile component and a term representing the correlation 
between U, and z.. For wind flow over. a rigid boundary, where the nature 
of the surface is indpendent of the flow parameters, the correlation 
term can be expected to go to zero and equation (3.8) reverts to the 
familiar logarithmic profile. We can interpret the first term on the 
right hand side of equation (3.8) as representing the mean velocity at 
any height required to support atmospheric dynamics, while the second 
term identifies profile modification due to external effects. In view 
of the data screening, the effects of errors in determining the profile 


descriptors should be expected to be random, rather than systematic 
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in the fashion shown in figure 3.5. Similarly, atmospheric stability 
variations are expected to be random in nature, and should not contri- 
bute to the correlation term. Richardson numbers were computed as 
instrumentation permitted, and indicated near-neutral stability condi- 
tions in the lower levels of the atmosphere over the ocean. Therefore 
the second term on the right hand side of equation (3.8) is taken to 
represent the wind at any height required to support air-sea interaction 
dynamics. For simplicity, we will write equation (3.8) as: 


U = Up +U, (3.9) 


Where U, is the profile component of the mean wind, and Uy is the boundary 
interaction component of the mean wind. Up is always a positive quantity 
however as shown in figure 3.5, U,, takes on both positive and negative 

1 


values, with a transition at about 2 meters/sec +. Given the mean para- 


meters of z, and Ux, the resulting momentum U, at any level z is all that 


Pp 
is required to effect a balance between the momentum created at some higher 
level and the momentum transferred or lost due to friction. The additional 


momentum, U., represents either the excess or deficit of momentum, where 


w? 
the positive sign indicates excess, available at level z to support the 
modification of the boundary characteristics. Figure 3.5 indicates that 
at wind speeds greater than 2 meters/ ee in the mean there is an 
excess of momentum in the atmosphere available for wave generation and 
maintenance. At wind speeds less than 2 neces cece in the mean the 


atmosphere does not appear to have adequate momentum to support frictional 


losses, and receives momentum from the ocean, presumably from the underlying 
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swell. This interpretation is consistent with and reinforces the view 

that there is a minimum wave generation wind speed at about 2 meters/ 
—] ~-1 ; : 

sec . Above 4 meters/sec the ratiovor U,/U varies between 0.1 and 

0.2. This ratio is consistent with the ratio of wave growth momentum 


to the total atmosphericmmentum determined by Stewart (1961). 


Conclusions. 
Charnock (1955) suggested z, should be proportional to Uy" /g 


resulting the wind profile formula: 


Ux Bz 
U(z) = = ee 
k Ux 
ee ) (3.10) 


8 is a constant of proportionality. Inspection of figure 3.2 clearly 
shows that z, does not vary as Ux*. Laboratory experiments by Wu (1968) 
resulted in similar findings, however Wu noted that the Charnock expres- 
Sion is appropriate at high wind speeds. This expression implicity as- 
sumes that a single surface roughness parameter will specify both the 
momentum balance requirements of the ocean and the atmosphere. There 

is no reason for this to be so. It is suggested that a more appropriate 
approach is the separation of the profile into two components, a compo- 
nent resulting from atmospheric dynamics and one resulting from wind 
wave interaction, similar to that done in equation (3.9). 

At low wind speeds, the vertical mean wind profile parameters are 
all functions of velocity, with the roughness length parameter being 
highly structured. Two of the three peaks in the roughness length graph 
can be explained by recourse to classical hydrodynamics. The peak at 


an ; : : 
4 meter/sec requires explanation. The rise and then the return to 
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a near zero value of the velocity defect in figure 3.5 at 4 mé ton ieccuan 
taken along with the roughness length plot suggests to this author that 
the peak at 4 meter eee may be associated with a transition from cap- 
illary waves to gravity waves as those waves having a dominant effect 
on the atmosphere. The evidence is,far too meager and the arguement 
far too tenuous, however, to make this assertion with any firm con- 


viction. 
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Chapter LV 
Small-Scale Atmospheric Processes Just Above the Sea Surface 


General 

Current theory describing the momentum transfer processes at the 
air-sea boundary require specific behavior of the atmosphere before 
momentum interchange can occur. A continuing controversy exists over 
the validity of the theory, based on the critical appraisal of the as- 
sumptions about the atmosphere directly above the water surface. Re- 
solution of the controversy lies in observation of the actual behavior 
of the atmosphere above the ocean surface. This chapter provides some 


observations in this regime. 


Background 


The wave generation mechanism proposed by Phillips (1957) relies 
on the advection of turbulent eddies over the sea surface to induce a 
resonant response at the ocean surface, resulting in the growth of wind- 
waves. Application of this theory requires a statement about how eddies 
are advected across the water surface. The theory uses the concepts of 
a dominant eddy convection velocity and an integral time scale of the 
turbulence in a convected frame. The theory describes the initial exci- 
tation of waves by turbulent air pressures which are not affected by 
wave motions. As soon as the waves start affecting the air motions, feed- 
back coupling occurs, which is described by Miles' (1957) theory. This 
theory predicts exponential ace caused by a vorticity interchange 
across a "critical level", defined as that level where U(z), the mean wind 


speed at height z, is equal to C, the phase speed of the growing wave. 





D6 


The rate of momentum transfer is proportional to 


: ve 
Maz (a*Tieeeicen eee (4.1) 


a is the density of air, while w is the vertical velocity. All quanti- 
ties in equation (4.1) are evaluated at the critical level. Application 
of the Miles theory requires an assumption about the constancy of the 
vertical mean wind profile and the magnitude of wt, The constancy of 
the mean wind profile appears valid as discussed in chapter III. The 
value of wie is usually taken to be equal to the mean square amplitude 
of the growing wave. This infers a coupling between wind and wave, 
where the motion of the wind above the wave is identically equal to that 
of the wave. 

Measurements of wave growth by Snyder and Cox (1966) and Barnett 
and Wilkerson (1967) indicate the theoretical wave growth is insuffi- 
cient to account for observed wave growth. Measurements by Seesholtz 
(1968), however, result in reasonable estimates of wave growth using 
Miles’ theory. In his results, Seesholtz demonstrated that wave-growth 
computations are quite sensitive to accurate profile measurements and 
that one must use mean wind speeds averaged over a short period of time. 

Stewart (1967) argues using wind spectra that waves do not induce 
fluctuations in the atmosphere of the order assumed by Miles of the we 
term in equation (4.1). Stewart's arguments lose force, however, due 


to the wide bandwidth of his spectral estimates and the lack of informa-~ 


tion about smoothing techniques employed. Hanning the spectrum, for 
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example, could spread the observed wave energy into adjacent bands, 
causing the observed low power at the wind-wave band. Weiler and Burling 
(1967) published spectra with narrower bandwidths. These spectra show 
considerable energy in the region of the typical wave frequencies; how- 
ever, sharp peaks were not observed. 

Some of the failure to observe wave energy in spectral density plots 
of atmospheric turbulence may be due to the method of representing the 
spectrum. To illustrate, let us examine a highly structured spectrum. 
Figure 4.1 shows two plots of an ocean surface wave spectrum, specifi- 
cally chosen because of the dominant peaks in the spectrum. Figure 4.1(a) 
is a linear plot of the wave spectrum, while figure 4.1(b) is a plot of 
f¢o(£) versus £n f£, as spectra usually are represented in atmospheric 
studies. Figure 4.1(a) indicates a dominant near-monotonic swell. Exa- 
mining figure 4.1(b), however, one would tend to conclude only that there 
is considerable energy in the region of typical swell frequencies. One 
does not clearly observe the sharp peak and the nature of the swell. Small 
detail is lost in the plot. In order to emphasize small features in the 
spectrum, all spectral density plots are presented herein with linear fre- 
quency scales. This is easily accomplished since only a narrow portion of 


the atmospheric frequency domain is presented. 


Taylor Hypothesis 


The Taylor, or "Frozen Turbulence" hypothesis relates the time and 
space scales of turbulence. This hypothesis postulates that if the mean 
wind speed were high enough, the turbulence would remain unchanged as 


it is being convected past a point, such that: 
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f=— (4.2) 


k denotes the wave number, U is the velocity of the mean air stream, and 
f is the frequency of the eddies measured at a point. In a shear flow, 
however, the Taylor hypothesis is expected to hold only for eddies of 


such size that: [see Lumley and Panofsky (1964 pg. 57)]: 


s LO (4.3) 





If we differentiate equation (3.1), and combine the result with equa- 
tions (4.2) and (4.3), the minimum frequency for which Taylor's hypo- 


thesis is expected to hold can be expressed as: 


Ux 


f >? 21KZ 





(4.4) 


where « is Von Karman's constant. 

Taylor's hypothesis has been repeatedly confirmed over land, as 
summarized by Monin (1967), resulting in its wide application over land 
and sea. The validity of the Taylor hypothesis in the frequency domain 
of wind waves has not been confirmed. In fact, the characteristic cats- 
eye pattern of the Miles wave-generation mechanism as described by Light- 
hill (1962) suggests that the atmospheric eddies would move coupled with 
the wind wave, rather than convected along by the mean wind. 

To measure the validity of Taylor's hypothesis, consider two wind 
measuring devices separated by a distance 2, and oriented such that they 


both lie on the same mean wind streamline. The average time Tt for 
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eddies to travel the distance © can be determined from the correlation: 
Ry? = <U(x,t), UGcra, t1)> (4.5) 


If tT is that delay time when R is a maximum, one can compute a correla- 
a y uu P 


tion velocity U,? 


eG a. 
Ue (4.6) 
a 


For Taylor's Hypothesis to be valid, Up =U. Table 4.1 lists the mean 


wind speed for two elevations measured by cup anemometers, along with 
U, as measured by hot film anemometers and determined from Equation 


R 


(4.6). Note that Up does not in general equal U. Further, there does 
not appear to be any pattern between the convected velocity and the ratio 
Ux/27Kz. 

For one day, the data allowed a more detailed analysis of the con- 


vected velocities. Equation (4.5) can be fourier transformed into the 


cross spectral density function, S(f), where: 


oO 


sf) =2 Ra (ena (4.7) 
eZ ope 
3 . (f) is a complex number and can be represented in complex polar 
ia 2 
form with the absolute value of Pups (f) and a phase angle $(f). This 


phase angle is related to sensor separation and the wavelength A of 


the ratio: 


o(£) _ 2m 
c i (4.8) 
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From Equation (4.2) the convected velocity can be expressed 


2T£L 
>(£) 





U,(£) = (4.9) 


In practice ¢(f) is measured only for values varying from 0 to 27. For 

A < S, however, the true phase angle has a value equal to ¢(f) + 2NT7. 

The criteria for determining the sensor separation, 5, did not foresee 
spectral computations. Consequently the value of S = 20 meters is too 
large and has resulted in the 2N ambiguity. This ambiguity has been ar- 
bitrarily resolved by selecting that value of UR nearest to the mean wind 
speed U, consistent with adjacent spectral values. In general, the 
correct value of Up was clearly evident. 

Figures 4.2 and 4.3 are two plots taken on the same day, where uN 
has been computed using equation (4.9) and plotted as a function of fre- 
quency. Both figures display the spectrum of the ocean wave field appro- 
priate for the measurement period; the measured mean wind speed for the 
appropriate heights; and the phase and group velocity for the underlying 
ocean waves. The values used to construct figure 4.2 were measured under 
calm, stable conditions. A heavy fog bank shrouded the measurement site 
and the ocean for a radius of at least 15 miles. Visibility at the site 
varied between 50-150 feet. U,/2mKz = .072 for 1 meter, .014 for a height 
of five meters. Figure 4.3, measured during the afternoon of the same 
day, was measured during strong wind-wave generating conditions. During 
this period, the sky was clear with unlimited visibility and estimated 
near-neutral stability conditions. U,/2mKkz = .11 for 1 meter, .021 for 


five meters. The spectral bandwidth in both figures is .05 Hz, the 
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standard error of estimate is .075 with 195 degrees of freedom. 

Examining the two cases, in the non-generating case (figure 4.2) 
those eddies corresponding to the frequency of the dominant waves ap- 
pear to travel faster than the same frequency eddies at 5 meters, while 
at high frequencies (f>.4) the eddies appear to be convected along at 
the speed of the mean wind. During the generating conditions (figure 
4.3) the picture has changed. The eddies at all frequencies except 
those at the dominant wave frequencies tend to move slower than the mean 
wind. 

This result suggests that the validity of Taylor's hypothesis is 
suspect in the near-ocean atmosphere that spectral domain encompassing 
frequencies at which ocean surface waves can be expected to occur. For 
all frequencies in figure 4.2 and for the swell frequencies in figure 
4.3, experimental information is inadequate to draw any inferences about 
possible phase lock, since there is no reason to expect the swell or the 
old wind-sea to be propagating in the same direction as the mean wind. 
The wind-sea in figure 4.3, however, was generated by the observed mean 
wind. In this case, we should expect the wind-waves to propagate in 
the same mean direction as the mean wind. Figure 4.3 suggests the exis- 
tence of phase lock for wind waves at both 1 and 5 meters. At frequen- 
cies greater than 0.4 Hz, there appears to be a tendency toward decoupling 
where some eddies move at a mean wind speed while some remain coupled 
with the wind waves. 

As discussed in chapter I, original plans called for the use of two 


remote buoys and a three-point space correlation field. The loss of one 
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of these buoys severely restricted a more definitive study of this pro- 


blem. 


Direct Reynolds Stress Measurements 

Application of the techniques discussed in chapter II resulted in 
the direct measurement of Reynolds stress for 5 cases. Unfortunately, 
a large amount of data was discarded since the data transmission reso- 
lution problem discussed in chapter II was not fully appreciated until 
late in the summer. Those runs of sufficient quality to process are 
listed in table 4.2. 

The tabulated drag coefficients agree with those of Weiler and 
Burling (1967), while the value of u'w' agrees with the square of the 
friction velocity, us, derived from the profiles. In case 5, the mea- 
sured value of u'w' appears high; however it is measured at a wind velo- 
city very close to that at which the boundary instability appears to 
occur as discussed in chapter II. Values of Ue may deviate from the 
observed u'w' in this regime due to the use of ten minute average mean 
wind speeds. Perhaps shorter averaging times are appropriate. [In all 


cases the normalized standard error of estimate is .10 (see appendix II). 


Wind-Wave Coupling and Capillary Waves 


During the course of the experiments, wind parameters were measured 
over a sea covered with an artificial sea slick. The measurements are 
compared with measurements made without a sea slick. These experiments 
provide detailed information about the near-surface wind field and pro- 
vide an indication of the importance of small waves to air-sea inter- 


action processes. 
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The sea slick was developed by Barger and Garret (1968) to provide 
intense capillary wave damping. Dr. Barger and Dr. Garret also developed 
the technique for dispensing the slick material, and worked directly with 
this investigator by providing the material and supervising its disper- 
sion over the sea surface. The experimental technique for applying the 
Slick and its gross results on both the wind field and the wind sea is 
described by Barger, et al. (1969). Briefly, measurements were made be- 
fore, during and after an artificial sea slick had been spread over the 
ocean surface. Marked modification of the wind profile roughness length 
parameter and the high frequency wave energy was noted. 

Using the techniques discussed in chapter II, the wind at a height 
of 1 meter above the mean water surface was resolved into a horizontal 
component, U, and a vertical eon e W. Each of these components, 
taken individually and together, were used to compute cross spectral 
density functions relating the variables to wave height as measured by a 
wave gauge. From these computations, the coherence between the spectral 
component of the wave field and the corresponding spectral component of 


the wind is computed. The coherence is defined: 


2 
B19 ‘*) 


wolle all (4.10) 
S11 (f) 5, (£) 


C(f) = 
where the spectral density function Si4 is defined in equation (4.7), and 
the subscripts 1 and 2 relate to any specified variables 1 and 2 respec- 
tively. 


Figure 4.4 displays the wave spectrum measured (a) well before the 
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slick arrived at the site, (b) just after the slick enveloped the site, 
and (c) just after the slick cleared the site. Note in particular the 
reduction of wave energy at high frequencies (greater than 0.3 Hz) as the 
slick persists. The increased wind-wave energy between the 1140 A.M. 
spectrum and the 1240 spectrum clearly identifies the condition as one 

in which wind waves were being generated. That portion of the spectrum 
above 0.3 Hz is the saturated portion of the wave spectrum, and should 

be constant from spectrum to spectrum. 

Figure 4.4 is constructed from a digitally computed spectrum. In 
order to eliminate uncertainty in the high frequency spectral estimates 
introduced by Nyquist folding, a spectrum was computed by analog tech- 
niques only for the high frequency components of the waves, using wave 
height data passed through a high pass filter and a waveform analyzer. 
Figure 4.5(a) represents the wave spectrum measured at the same time 
the 1140-1210 spectrum in figure 4.4 was measured. Figure 4.5(b) was 
measured just before the slick cleared the measurement site, while figure 
4.5(c) shows the wave energy well after the slick had cleared the measure- 
ment site. The difference in high frequency wave energy from spectrum 
to spectrum is clearly evident. 

We will compare the wind field over the ocean during the time the 
sea slick was in the vicinity of the measurement site with the wind field 
observed after the slick had passed. This choice is based on the simi- 
larity of the ocean surface wave spectrum for each period as shown in 
figure 4.4. The wave spectrum measured before the slick arrived at the 


site differs enough from the other two spectra to eliminate this run for 
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direct comparison of the data. 


The cross spectral density and coherence was computed for the fol- 


lowing variables measured during and after the passage of the artificial 


sea slick: 


a. 


b. 


Horizontal velocity, U, versus wave height, n. 
Vertical velcoity, W, versus wave height, n. 
Horizontal velocity, U, versus vertical velocity, W. (Reynolds 


stress). 


Plots of this data are displayed in figures 4.6 through 4.8. For all 


plots, the spectral bandwidth is .025 Hz, and the standard error of es- 


timate is 0.10 with 100 degrees of freedom. These plots indicate: 


a. 


The coherence between the swell and the wind field is modified 
only slightly, if at all, by the presence of a sea slick. With- 
out the slick, however, the horizontal velocity maximum occurs 
over the wave trough. With the slick present, the horizontal 
velocity maximum shifts to over the wave crest. 

The coherence between the wind sea and the wind field directly 
above the water surface is strongly influenced by the presence 
of the slick. With the slick absent, there is a moderate coupling 
of the wind and the sea. When the slick is present, the wind 
field appears to lose all knowledge of the motions of the wind- 
sea. This can be seen in the U and W versus n coherence plots 
and also in the U versus W phase angles. The phase angles of 
this Reynolds stress spectrum are completely unaffected by the 


underlying sea during the presence of a slick. 
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c. The horizontal velocity (figure 4.6) is a maximum over the wave 
trough for both the wind sea and the swell. This result was 
first observed by Seesholtz (1968). At that time some had argued 
that such might be true for either the wind sea or the swell; 
however, figure 4.3 shows this phase relation to be true for 
both. 

On the basis that the primary effect of a sea slick is the damping 

of capillary waves, the experimental evidence suggests that the coupling 
of the long wavelength wind-waves to the atmosphere is controlled in some 
manner by the capillary or short wavelength gravity waves. This would 
suggest that theoretical models for wave generation need to address both 
the long wavelengths generating wave and its relation with the short 
wavelengths wave. Similarly, it would appear that laboratory investiga- 
tions may be of limited value in wind-wave generation studies unless 

both the long wavelength gravity wave and the short wavelength capillary 
wave can be simultaneously scaled and represented in a laboratory environ- 


ment. 
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Chapter V_ 


Sumnary of Results and Recommendations 


The Direct Measurement of Reynolds Stress 


The techniques outlined in this report appear to provide valid 
results for the direct measurement of Reynolds stress. Unfortunately, 

a large portion of the summer effort was expended perfecting the techni- 
que, so the results were not as copious as desired. The data presented 

in chapter IV represents measurements made under near ideal equipment 
conditions. That is, all instrument gains were set at near optimum levels 
for data transmission, while system noise was low. The measurement system 
appears to be capable of surviving in the environment, if one can insure 
faithful data recording and conversion. 

Availability of digital processing or sophisticated analog data pro- 
cessing is mandatory for processing the data from the Reynolds stress 
sensor. The individual probe signals should be recorded either analog 
or digitally without linearizaton and then linearized by the terminal 
processor. If digital processing is used, a minimum of 12 bits of signi- 
ficance is considered necessary to avoid loss of data significance due 
to quantitizing errors. Significance was achieved with a 10-bit digi- 
tizer only through careful data selection and gain adjustment. 

It is recommended that continued direct measurements of Reynolds 
stress be made, but only after facilities are available to provide digi- 
tal data with 12 bits or more of significance for processing. The pre- 
ferred alternative is digital recording in the field. This step by passes 
intermediate analog recording and processing, eliminating the attendant 


opportunity for data degradation. 
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f 
The DISA anemometer systems used for measuring Reynolds stress 


proved to be ideally suited to the task. The weak point in the system 
is the step resistance changes of eremereneee This probe behavior is 
unsatisfactory and needs correction. It would appear that the probe 
behavior is more a result of manufacturing quality control than design. 
A probe purchased in 1967 performed for two years without trouble, 
while all of the probes purchased in 1968 suffered unacceptable resis- 


tance variability. 


Vertical Mean Wind Profiles 
Details of the vertical mean wind profile parameters are beginning 
to emerge. In the wind regime studied, the wind speed and the friction 
velocity are linearly related in the mean, while the roughness length 
parameter is highly structured. Evidence suggests a minimum 10 meter 
wind speed for wave generation of 2 meneea eeu The roughness length 
and friction velocity plots indicate the existence of boundary instability 
at a ten meter wind speed of 4 and 8 meeen/ Soome The instability at 
8 pesicec can be explained by the Kelvin-Helmholtz instability. 
The one at 4 meters/sec - needs explanation. Finally, the mean profiles 
were shown to have a component which contributes to boundary modification. 
It is important to consider the profile results in terms of an air- 
sea momentum budget. That momentum lost by the atmosphere can be partitioned 
into a portion lost to friction, a portion lost to the maintenance of the 
existing wave field, and a portion required to generate new waves. Gther 
effects, such as drift current generation, can be lumped together with 


the friction term. Use of the profile parameter plots as models of the 
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atmesphere in some wave-generation prediction scheme may be ill-advised 
since the plots do not provide for the partition of momentum between 
that required to maintain the wind sea and that provided for further 
growth of wind waves. Failure to consider the momentum partition may 
have further intorduced a systematic bias into the plots. At intermediate 
wind speeds one might expect generating conditions to be observed more 
frequently while at high wind speeds the momentum required for maintaining 
the existing wave field may become the dominant effect. Therefore the 
plots in chapter III are presented as representing a "typical" mean 
state of the vertical mean wind profile over the sea. 

Clearly, the next step in the study of the vertical mean wind 
profile over the sea is the consideration of momentum partition. It 
is recommended that observed profiles be catalogued and sub-divided into 
three sets; one set containing profiles observed during rising seas, one 
set containing profiles observed during falling seas, and finally one set 
of profiles observed during fully arisen seas. Each set can be analyzed 
using the techniques followed in chapter III. The resulting profile 
parameters will represent mean prufiles appropriate for use in modeling 
an air-sea relationship. 

In order to provide for the fine grained analysis of mean wind 
profiles over the ocean, the profile data set must be increased by 
at least a threefold amount. It is recommended that efforts continue 


to collect mean wind profiles over the ocean from rigid platforms. 


Small scale wind motions above the sea. 


The Taylor Hypothesis does not appear to be a valid approximation 


for treating turbulence directly above the ocean surface in the spectral 
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domain of ocean surface waves. 

Sea slicks whose effects on ocean surface waves is the intense 
damping of capillary wave systems are shown to have strong effects on 
the coupling of the near-ocean wind field with the long wave length 
ocean wind waves. This result suggests that the capillary waves play 
an important role in the generation of wind waves and the dynamics of 
the atmospheric boundary layer above the ocean surface. 

The horizontal velocity maximum occurs over the wave trough, both 
for the swell and the components of the wind sea. 

These results demonstrate the importance of observing small scale 
atmospheric processes near the ocean surface. The use of artificial 
sea slicks provides the field investigator with a powerful research 
tool. Their use constitutes one of the few measures a field researcher 
can employ to control his environment. It is strongly recommended that 
this technique be exploited to its fullest to enhance our understanding 
of the small scale processes at the air-sea boundary. 

The apparent importance of capillary waves and/or short gravity 
Waves to air-sea dynamics is a matter requiring additional investigation, 
both theoretical and experimental. Before realistic models of the air- 
sea boundary processes can be constructed, the role of these waves in 


the total dynamics must be identified. 





ion 


Acknowledgements 


Professor Erik L. Mollo-Christensen at the Massachusetts Institute of 
Technology has been my principal source of advice, support and instruc- 
tion during the past several years. His encouragement and patience is 
most sincerely appreciated. 

Mr. Ortwin von Zweck worked closely with me both in the field and 
during the analysis phase of the project. Our many long discussions 
and hours working together on common problems has been a key to the 
success of this study. 

Mr. Kenneth A. Morey and Mr. James H. Peers worked on the technical 
end of this study since its inception. Their innovative skill, techni- 
cal competence, and willingness to devote many long hours to the support 
of this work contributed to the success of the project. 

There are many who have contributed each in his special way to this 
effort. Mr. Anthony Cella worked many hours with me on Lewis Wharf assem- 
bling the prototype satellite buoy system during the early stages of 
this project. Mr. Edward Bean's skill in the machine shop overcame 
many of the design problems in the gaeetuftes buoy. Mr. Louis LeTourneau 
operated the research boat during the field portion of the work. 


Miss Louise Harris and Miss Diane Lippincott typed the manuscript. 





ee 


Appendix I 


Wind Profile Data 


Elementary statistical techniques have been applied to the raw 
wind profile data to develop this table of data. The usual form for 


the logarithmic wind profile can be written: 


K 


Te U(z) + &n(zo) : (APD) 


Qn(z) = 
where z is the height above the boundary, U(z) is the mean wind at 
height z,Ux is the "friction velocity", where Uy? = - U'W'; K is von 
Marman's constant and is equal to approximately 0.42 and z, is the 
roughness length. Equation (AI.1) is linear in 2£n(z), therefore the 
logarithmic wind profile can be fitted to the raw data using linear 
regression techniques (see chapter 3). Such a best fit establishes 
the form of the profile, the roughness length, z,., and the friction 
velocity, Ux. gee a correlation coefficient can be computed 
providing a measure of how accurately the logarithmic profile repre- 
sents the raw wind data. 

| After the form of the profile had been established, the wind data 
was printed for the "standard'heights as shown in the table. When raw 
wind data was available for the standard height, it was printed, other- 
wise a wind velocity computed from equation (Al.1) was printed. The 
standard observation heights for all data during the summer of 1968 was 
1 meter, 2 meters, 3 meters, 5 Meee een and 8.2 meters. The observation 
heights for Seesholtz' data are as published by Seesholtz (1968). 

The drag coefficient is defined as: 


ee * )? (AI. 2) 


= ell =p 


ad 
eee, o-e 


is 
aa Sy a 
- 
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80. 


This was computed using the value of friction velocity derived from 
the profile for 2 meters and 10 meters. 


The gradient Richardson number is defined as: 


30 
ee Oe 
Ri = 
ay? (AL. 3) 


For purposes of computation, equation (AI.3) was rewritten: 


on Saas At AU 2 
Ri = (=) ce rD (1/¢ eo (AT.4) 


Where T was the mean temperature at 3 meters, AT/Az the temperature 


gradient between 1 and 5 meters, I, the dry adiabatic lapse rate, 


d 
and AU/Az the wind gradient between 1 and 5 meters. g is the accel- 
eration due to gravity. 

Each profile is identified by the month, day and run number. For 
example, profile 08-14-203 is the third profile taken during run 2 
on 14 Augst 1968. The profiles carrying an "S" in front of the profile 
identity were compiled from Seesholtz' data. 

Further details on the measurement and computation of the pertinent 


profile information is contained in the appropriate sections of this 


report. 
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Appendix II 
Digital Data Processing 


Initial plans for the collection and analysis of field data called 
for the use of analog data processing, using procedures similar to those 
described by Seesholtz (1968). It became clear, however, that the type 
of data processing required to adequately investigate the physical phen- 
omena demanded the use of digital computation. The principal limitations 
of the existing analog system were: 

1. The inability to compute cospectra and quadrature spectra. 

2. The distortion of spectral estimates caused by linear trends 
in the data. This problem was most pronounced in the wind data, where 
Strong linear trends were often present. 

3. The inability to compute probability density and joint prob- 
ability density distributions. 

4. The restricted precomputation and preprocessing of data before 
Statistical measures were computed. 

The shift to digital processing introduced a new set of problems 
in the overall data handling. The major digital processing handicaps were: 

1. The need to develop the processing programs from scratch. 

2. Analog-digital conversion facilities with a maximum 10-bit 
resolution. A minimum of 12 bits are required for significance with 
some of the data (see Chapter II). 

The programs were completely written and tested out in a period 
of two months. This was accomplished byusing the most direct analytic 


methods possible in all computations. A more refined program could have 
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been written, however the time required to master the subtleties the 
computer or advanced computation techniques was not considered warranted. 

The analog to digital conversion facility is maintained by the 
Department of Mechanical Engineering at the Massachusetts Institute of 
Technology. The facility consists of a Raytheon 16 channel Analog to 
digital converter connected to an IBM 1130 computer. The converter 
will accept signals with a maximum signal level of 110 volts, and will 
sample at any preset sample rate up to about 6000 samples per second. 
The final output is standard IBM punched cards formatted with 56 data 
points per card. 

Analog data tapes were prepared for digitizing by prefixing each 
data tape with a leader containing a square wave calibration signal. 
The tapes, along with their leaders, were then digitized at 10 times 
the recorded tape speed. The calibration signal at the beginning of 
each tape allowed optimizing the gains of each channel amplifier into 
the digitizer to fully use the 110 volt digitizing window. After IBM 
cards had been punched for each data channel, the first several hundred 
points, which included the calibration signal, were printed out. The 
digitized value of the calibration signal was then used to determine a 
scale factor for the data to be used in later computation. 

The general procedure for digital processing followed those recom- 
mended by Bendat and Piersol (1966). The theoretical maximum sampling 


period is: 


rae (AIL.1) 








oi 


Where E. is the maximum significant frequency in the record. For wind 
studies over the ocean, there is little energy of contern above 0.5Hz, 
however on occasion there may be reason to study up to 1.Q0Hz. For 
accurate correlation function measurements, however, Bendat and Piersol 
recommend: 


at = aE | aan 


On this basis, most of the records were digitized with a digitizing 
period of 0.2 seconds, while some were digitized with a period of 
0.4 seconds. 

While correlation functions require data with high sampling rates, 


the resolution bandwidth for power spectra is given by: 


e ~ Atm (AII.3) 
where m is the number of correlation lags. Similarly the normalized 


standard error, €, is specified by: 


m™m 
S -/2 (AII.4) 


Where N is the sample size. The computation cost is directly proportional 
to m, while the error increases with increasing m. It is therefore ad- 
vantageous to specify as small a value for m as possible. From equation 
(AII.3), however, the spectral resolution decreases with decreasing m, 

and increases with both increasing m and At. Consequently for spectra, 

as large as possible At is desireable. The program has been set 


up, therefore, to provide the basic At for computations requiring 
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high sample rates, however the sample period At is doubled by discarding 
every other point f-r spectra computations. 

After the analog data was digitized and the cards prepared for 
processing, the following digital computations were performed in the 
order listed: 

1. Input. Data, control information, and supplementary information 
was read into memory. The program could handle four simultaneous channels, 
each with 9600 data points. 

2. Scaling. Each channel of data was multiplied by a scaling con- 
stant provided in the input to restore the digitized data to physically 
meaningful values. 

3. Special processing. Special processing included manipulation 
of data as required for this study and other program users. Special 
processing included the resolution of hotwire data as described in 
Chapter II, and the building of data array consisting of the product 
U(t)W(t) for correlation of Reynolds stress with wave height. 

4, Computation and removal of mean and trend. The mean was computed 
using the formula: 

x =e Jee (AII.5) 
es 
For computing the trend, cognizance is taken of the fact that the indexing 
variable, i, is a function of time since each adjacent storage location 
is filled with data sampled at even intervals of At. Consequently 


one can compute the linear regression formula where the value of the 
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variable is correlated with its locetion in the array. Therefore, the 


predicted value of x at any point (denoted as x')is from (Hoel (1954, 


N 


pel27)): ; 
: ae N 
i=, “ii - N x (5) N i 
x' = ibs Lae Ci =) ) +x (AII.6) 
> ih () 
i=l 


Simplifying, and noting that we want to subtact x' from each value of 


; the corrected value of the data in each array, x, becomes: 


x, le: X, ~ X is 


i = 5 x 
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The quantity in the brackets in equation (AII.7) is the slope of 
the linear trend line. This quantity is printed out with the other 
Statistical measures to allow the user to assess the magnitude and effect 
of the trend. This technique of trend correction eye a clear 
Spectrum of spectral components of frequencies equal to or higher than 
fy; = 1/N(At). Estimates of the power spectral density between 0 and fy, 
will be in error unless there isn't a trend. The low frequency com- 
ponents will appear as a higher frequency. This can best be graphically 
illustrated. Figure AII.1 shows an arbitrary frequency component of 
f<f;, with period T. The trend line, AA', is drawn as the best fit of 
this component, resulting in an apparent period T'. The period T' is 
less than T, showing that the low frequency component is represented 
as a higher frequency. Consequently, while the technique of linear 


trend removal is simple and direct, the first data point of the raw 
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Figure AII.1 


Frequency shift of low frequency components 
caused by linear trend line AA'. Frequency 
component with period T is represented as a 
frequency with period T'. 
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_ Spectrum or the first two data points of the smoothed spectrum will have 
incorrect power spectral density estimates unless the linear trend line 

has a slope near zero. In practice, since the maximum array capicity 

is 9600 data points, a slope of order 10s or smaller can be interpreted 


as a zero slope. 


5. Computation of higher order moments. The second through the 


fourth moments are computed and listed for the data in each channel 


using the following formulas: 


aay 1 2 2 
variance X“ = Wo 2 x. ] 
i=l] a 
N 
‘ieee eee ei 
5 Nae i 
i=] 
Li Summee a7 
Kurtosis M, = [ N z xX. ]) (X*) 


i=l (AIL.8) 

6. Probability Density. A frequency histogram or a joint frequency 
histogram for any variable or any pair of variables can be constructed 
for up to a maximum of 40 dates intervals for each variable. The 
technique is an elementary sorting procedure. 

7. Correlation function. Following Bendat and Piersol (1966), 
the correlation function can be computed for M lags: 


i=N-} 


1: 
ey) REY _ a iat j= 0,1, 2, .-.M  (CAII.9) 


Where the subscript j denotes the jth lag, corresponding to a time lag 


of jAt. For jAt = tT, equation (AII.9) may also be represented as 
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BoD = <x(t)y(ttT)>. If Y =X then equation (AII.9) represents an 
autocorrelation where: 


Sw = Sa (ATT S10) 


Consequently when x = y only the function Se is computed. For a 
cross-correlation, however, Rey 6) and ae is computed, where 
ae” = Bt 

8. Power spectral density. The final computation is the fourier 
transformation of the correlation function into the cospectral density 
function and the quadrature spectral density function and then repre- 


senting these as a power function and phase angle. The correlation 


function can be represented as an odd and an even part: 


Aj + (RG) + R bs (j)] 
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For an autocorrelation, Aj = Seg I and Bj = 0, however for a cross 
correlation, both Aj and Bj may be different from zero. If Sg Be 
is the co-spectral density and Quy SF) is the quadrature spectral 


density, then the spectral density estimates are: 


Moai 


GC. (£) = 2At{[A, + 2 3% Aj cos ( a1 ) + Ay cos (]{M£)] 
jal Cc £ 
Cc 
M-1 
Q(£) = 2ae[2 2 By sin AL) +B, sin re) 
a= c ConeCAgolieslio)) 


where: 
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Smooth estimates are calculated by smoothing the spectral density 
estimates by the Hanning method: 


Cy = 0.25 Ca 7025 Cy tc Ohezs) Chad 


OMe 02 5a Ose (Orme) Sunes 


Ce )0 5 io 025 Cy 


Onn 0. Dae tt Oe Q) 
Cc. = 085 © + 0.5 € 
m m-1 


Onan a0. Omae (AIT. 13) 


Finally the spectral density estimates can be represented as a power 


and phase angle: 


7 2 2 
Py (f) = (Cy + Q?) 
O., (£) = tan “4 (2x) : (AII.14) 
k 


These are listed and plotted for analysis. 
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Appendix III 


Hotwire Linearizer Design 


General Discussion 

Effective field use of hotwire anemometers for the study of meteor- 
ological turbulence is complicated by the inherently non-linear output 
of the anemometer. Of the several solutions to this problem, the simplest 
and most direct is to linearize the hotwire anemometer signal in the field 
before recording. This approach, from a data processing viewpoint, facil- 
itates system calibration, error tracing, and bookkeeping while recording 
and processing the data. 

While the advantages of hotwire linearization are usually self-evident, 
economic dictates usually result in adopting alternate procedures. The 
linearizer discussed herein is one approach to sea idinre direct hotwire 
linearization while circumventing the cost problem. The linearizer has 
limited performance, and is unsuited for high frequency measurements. 

High frequencies are seldom of significance in the field study of meteor- 
ological turbulence, however, making the instrument suited for this appli- 


cation. 


Design 
The response of a hotwire system can be represented by: 

E” - Eo” = RCYU Gam 
Where E is the output voltage of the anemometer, E, is the offset volt- 
age, C the response characteristic of the anemometer, R the operating 
resistance of the hotwire, and U the wind velocity. The output voltage 
E can be written: 

E =) Be +748 (A3.2) 


where AE is the difference voltage between the offset voltage and the 
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output voltage. Substituting equation (A3.2) into (A3.1), the neglecting 


terms of order AE*or greater, the transfer function of the linearizer 
becomes: 
2 


2 
SE eer: (A3.3) 


CR 





US 
The terms enclosed by the parenthesis are constant for any given anemometer 
system, and need only be set into the linearizer during the initial cali-~ 
bration of the linearizer. 
In order to express equation (A3.3) in a form suitable for analog 
computation, the products are expressed as a sum of logarithms: 
Log K+ 2 log AE = log U (A3.4) 


where: 


2Eo ye 
CR 


oe 

Equation (A3.4) can be represented in terms of amplifier building blocks 
as shown in figure 1. Using standard non-linear amplifier circuits, 
the circuit as shown in figure 2 was developed. 

The construction of the linearizer is straight forward, using 
standard components. 

The functional relationshipsincorporated into the linearizer crucially 
depend upon the accuracy of the resistors used in the construction of 
the set. Consequently all resistors used in the linearizer should be 
temperature compensated precision metal film resistors. 

The frequency response of the linearizer is controlled by capacitors 
Cl and C2. The frequency response of the linearizer can be extended by 
replacing these capacitors by smaller valued units, say .001 microfarads. 
Care should be excercised here, however, since excessive reduction of 


the values of Cl and C2 will degrade the stability of the linearizer. 
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Within limits, one can tailor the high-pass characteristics of the lin- 
earizer through judicious selection of Cl and C2 to best fit experimental 
requirements. 

While equation (AIII.4) shows U to be a function of AE*, the higher 
order terms become important at high wind velocities. This along with 
the fact that some hotwires have response functions which deviate from 
the ideal King's law response, results in a need for a variable exponent. 
Consequently a multiple exponent selection has been incorporated into 
the design, providing a choice of exponents of 2, 2.5, 3, and 4. The 
response curve of the linearizer appropriate to the test system used in 
conjunction with the hotwire set is shown in figure 3. The frequency 
response of the system is shown in figure 4. 

Since a number of these linearizers were required, printed circuit 
boards were drilled and etched. This construction technique provides 
mechanical rigidity and compactness, in addition to the ease of construc- 
ting multiple units. Careful layout of the circuit and mechanical rigidity 
of the unit is necessary to guarantee circuit stability, particularly if 


the unit is destined for exposed field use. 


Operation. 


The output of the linearizer is a negative voltage varying over 
a range of approximately zero to five volts. The relationship between 
the output voltage and the input voltage is a function of the exponent 
and the operating range chosen, and must be determined by set calibration. 
To initially set up the linearizer for use, the following steps 


are suggested: 
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1. Determine E, for the hotwire set through wind tunnel calibration. 


The standard technique is to plot the hotwire set response on a graph 
with Re as the abscissa and VU as the ordinate. The best fitting straight 
line is drawn through the plotted points, and extended until it crosses 
the voltage axis. The vaiue of Bee is then read off at the intersection 
of the voltage axis and the response line. 

2. Adjust the linearizer for E,. Impress a voltage across the 
input of the linearizer equal to E,. Connect a high impedence voltmeter 
to test point A and ground. Adjust R1 until the voltmeter reads about 
0.2 volts. Theoretically, one should adjust this voltage to read zero 
volts, however the logarithmic amplifier tends to depart from a logarithmic 
function for voltages near zero, and a small “keep-alive" voltage is 
kept across the amplifier to improve the response characteristics. 

3. Adjust the linearizer scale. Set the exponent to the desired 
range. If the exponent value is not known, try 2.5 as a first approximation. 
Impress a voltage across the input equivalent to the hotwire voltage 
output corresponding to the maximum expected wind speed. Adjust R2 for 
a linearizer voltage output of 5 volts or any other convenient voltage 
depending upon the range of wind speeds of interest. 

4. Check and @librate linearizer. Run a calibration curve on the 
linearizer to determine the suitability of the exponent and range setting 
Readjust as necessary. 


9. Connect the linearizer to the hotwire anemometer and operate. 


Critique. 
The linearizer has obvious drawbacksin both circuit performance 


and in the ease of setting up the instrument for use. This tradeoff is 


necessary to design an inexpensive unit. Ideas come to mind where the 





106, 


ease of operation and the performance could be improved with incremental 
cost increases, however successive increments of improvement results in 
a deviation from the basic design objective to providing an operational 
linearizer at the cheapest possible cost. On the positive side the 
linearizer herein discussed runs at least ten times cheaper than most 


commerical models. 
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